NPS69-88-005 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




PERFORMANCE COMPARISONS FOR TWO VERSIONS 

/■ 

OF THE STANIFORTH-MITCHELL 3AR0TR0PIC 
NUMERICAL WEATHER PREDICTION CODE 

R. E. NEWTON and A. L. SCHOENSTADT 
" July 1988 



Interim Report for Period 
October 1987 - December 1987 



Approved for Public Release; Distribution Unlimited 
Prepared for: 

m 1 ■'nvi ronmental Prediction Research Facility 

/ , California 93943 



FedDocs 
D 208.14/2 
NPS-69-88-005 



DUDLEY «r>.TOX LIBRARY 

NA V Ai -).Sm :R. )I. ATE SCHOOL 

MONTEREY, CALIFORNIA S3943-5002 

f ; 17 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral R. C. Austin Harrison Shull 

Superintendent Provost 



This report was prepared in conjunction with research 
conducted for the Naval Environmental Prediction Research 
Facility and funded by the Naval Postgraduate School. 



Reproduction of all or part of this report is authorized. 



This report was prepared by: 



CURITY Cl ASMMCAHON 01 iMlSPAGt 



DUDLEY KNOX LfSRARY 
hWVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



REPORT DOCUMENTATION PAGE 



la REPORT SE^URI 

Imciassi fied 



TY CLASSIFICATION 



lb RESTRICTIVE MARKINGS 



2a SECURITY CLASSIFICATION AUTHORITY 



2b DECLASSIFICATION / DOWNGRADING SCHEDULE 



3 DISTRIBUTION /AVAILABILITY OF REPORT 

Approved for public release; distribution 
is unlimited. 



4 PERFORMING ORGANIZATION REPORT NUMBERS(S) 

NPS69-88-005 



S MONITORING ORGANIZATION REPORT NUMBE R<S) 



6a. NAME OF PERFORMING ORGANIZATION 
ORGANIZATION 

Naval Postgraduate School 



6b. OFFICE SYMBOL 

(If rfppitCdbJe) 

69 



7a NAME OF MONITORING ORGANIZATION 

Naval Environmental Prediction Research 
Facil ity 



6c ADDRESS (City. State, and ZIP Code) 

Monterey, CA 93943 



7b ADDRESS (City. State, and ZIP Code) 

Monterey, CA 93943 



8, NAME OF FUNDING/SPONSORING 
ORGANIZATION 

Naval Postgraduate School 



Bb OFFICE SYMBOL 

(if Applicable) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 

O&MN, Direct Funding 



8c ADDRESS (City. State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 





PH OGham 


PHOJI CT 


I ASX 


WOHk UNIT 


Monterey, CA 93943 


(It MINI NO 


*0 


NO 


ACOSSlON NO 



1 1 TITLE (include Security Classification) 

PERFORMANCE COMPARISONS FOR TWO VERSIONS OF THE STANIFORTH-MITCHELL BAROTROPIC 



12 PERSONAL AUTHOR(S) 








R. E. Newton, A. L. 


Schoenstadt 






13* TYPE OF REPORT 


13b TIME COVERED 


14 DATE OF RE PORT (tc^r Month. iMy) 


15 PAGE COUNT 


Interim 


from .10/1/87012/31/87 


26 July 1988 





16 SUPPLEMENTARY NOTATION 



17 



COSATI CODES 



FIELD 


GROUP 


SUB-GROUP 















1 B SUBJECT TE RMS (Connnoe on reverse tl netevwfy tacntily by block number) 

Finite element; numerical weather prediction 



9 ABSTRACT (Continue on reverse if rve<e%wr> and tdenltfy by block number) 

In a previous report R. E. Newton proposed \two improvements to the Staniforth-Mitchel 1 
jarotropic numerical weather prediction code. Reported here are results of perform- 
ince comparisons between the original code and an amended version which incorporates 
the two improvements: a new solution routine for the ei genprobl ems and a direct solu- 
;ion of the Helmholtz equation. It is found that the proposed improvements do reduce 
:he computation time, but that the new eigenproblem algorithm results in excessive 
•ound-off error when single-precision arithmetic is used. 



I DISTRIBUTION / AVAILABILITY OF ABSTRACT 
UNCLASSIFIED /UNLIMITED □ SAME AS RPT QDTICUSERS 



21 ABSTRACT SECURITY CLASSIFICATION 

Unci assi fied 



a. NAME OF RESPONSIBLE INDIVIDUAL 

. E. Newton 



22b TELEPHONE (irtciuoe Ace# Code) 

408 646 3424 



22c OFFICE SYMBOL 

69 



FORM 1473, 84 mar 



8 J APR edition may be used until exhausted 
All other editions are obsolete 



cy*n » uasw it a i 19*1 o* ih?> ** W - . 

# IJ S PtmiifVT OMm p lVbb *Ot> >4 




■'’OCvop a t-vyt a 



^ A Err ATE SCHOOL 

1 cv^cm^ t 



PERFORMANCE COMPARISONS FOR TWO VERSIONS OF THE 



' lJX ^3943-500: 



STANIF0RTH-M1TCHELL BAR0TR0P1C 



NUMERICAL WEATHER PREDICTION CODE 



I ntroducti on 

Reference 1 proposed improvements to the Stani f orth-Mi tche 1 1 
(Ref. 2) barotropic finite element code. It was anticipated that 
these improvements would effect significant reductions in compu- 
tation time. Reported here are comparisons of calculated results 
and times required to perform key sets of calculations using the 
"original” code with corresponding times for an "amended" code 
which incorporates the proposed changes. The test problem uses a 
periodic East-West boundary condition, walls on the North and 
South boundaries, and a 12x13 grid. Timing data were obtained 
using an IBM 4381-MI processor operating in the batch mode. (It 
was found that the system utility routines SETIME and GETIME do 
not give useful results when used in the time-sharing mode.) 

Timing Results 

E i genpr ob 1 em . The first comparison is for the determination of 
eigenvalues and eigenvectors for the x-direction "stiffness" and 
"mass" matrices. In the original program the relevant routines 
are MTRXC and EIGEN2, both called in EBVSET. In the amended 
version the routines are SETABX, SETD2N, and PEREIG, all called 
by SETUP. Within PEREIG there are 3 calls to EIGVCP and 2 calls 
to RAYLYP. The time for the original version is 104 milliseconds 
and for the amended version is 34 milliseconds. Because the eig- 
enproblem is solved 4 times in the original program and only 3 
times in the amended version, the overall comparison is between 
415 milliseconds and 101 milliseconds. 
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T i me Step. For each time step the Helmholtz equation is solved 
once. In the original version an iterative scheme devised by 
Concus and Golub is used, whereas the amended version uses a 
direct solution. When only a single iteration is employed in the 
original version, the difference in solution times is entirely a 
result of the extensive rearrangement of the right-hand side of 
the Helmholtz equation that occurs in the Concus and Golub 
scheme. This rearrangement is effected in SOLHEL. It is believ- 
ed that the most useful comparison is between the times required 
for a single execution of the routine TSTEP. For the original 
version this time is 372 milliseconds and for the amended version 
it is 331 milliseconds - a saving of 11%. 

Remarks 

In producing a version of the barotropic code incorporating 
the changes proposed in Ref. 1 it was discovered that calculated 
results diverged from those of the original version. It was not 
immediately obvious which results were better. The likely cause 
seemed to be differences between the two solutions to the eigen- 
problem. A comparison was made by constructing the product of 
the back-transformation matrix with the forward transformation 
matrix for each version. Since this product should yield the 
identity matrix, measures of the quality of the transforms were 
obtained by finding the standard deviation (from unity) for the 
elements of the principal diagonal and the standard deviation 
(from zero) for the off-diagonal elements. For the diagonal ele- 
ments the magnitude for the original version was 3.5E-06 and for 
the amended version it was 5E-06. For the off-diagonal elements 
the disparity was greater: 5E-07 for the original and 25E-07 for 
the amended version. 
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Since the e i gens o 1 u t i on of the amended version is iterative- 
ly improved by successive calls to RAYLYP and E1GVCP, efforts 
were made to improve results by additional calls to this pair of 
subroutines. Although this did result in addition of one cycle 
beyond the single cycle recommended in Ref. 1, no further im- 
provement seems to be possible using single-precision arithmetic. 

In order to resolve doubts concerning the e i genso 1 u t i on 
scheme of the amended version, a separate double-precision pro- 
gram was written to compare the eigenvectors and eigenvalues with 
exact results. Using 3 cycles of improvement, both eigenvectors 
and eigenvalues were accurate to 15 decimal digits. 

Cone 1 us ions 

In Ref. 1 attention was directed to the marginal accuracy 
of single-precision arithmetic for this program. The present 
study seems to have encountered further evidence in support of 
this observation. Although the amended version shows some reduc- 
tion in computation time, its adoption cannot be recommended 
without first verifying, by use of double-precision arithmetic, 
that it does indeed give results compatible with those of the 
original version. 

While the present study was proceeding, a new possibility 
for major improvement in computational efficiency came to atten- 
tion. In Ref. 3 Temperton and Staniforth describe a new time in- 
tegration algorithm which they call semi -Lagrangi an, semi-implic- 
it. They show that it allows much longer time steps than the 
present scheme while still giving equal accuracy. Serious con- 
sideration should be given to incorporating the new algorithm in 
the barotropic program. 
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